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ABSTRACT 



Context. Most exoplanet imagers consist of ground-based adaptive optics coronagraphic cameras which are currently limited in con- 
trast, sensitivity and astrometric precision, but advantageously observe in the near-infrared window (1 - 5ytim). Because of these 
practical limitations, our current observational aim at detecting and characterizing planets puts heavy constraints on target selection, 
observing strategies, data reduction, and follow-up. Most surveys so far have thus targeted young systems (1 - 100 Myr) to catch the 
putative remnant thermal radiation of giant planets, which peaks in the near-infrared. They also favor systems in the solar neighbor- 
hood (d < 80 pc), which eases angular resolution requirements but also ensures a good knowledge of the distance and proper motion, 
which are critical to secure the planet status, and enable subsequent characterization. 

Aims. Because of their youth, it is very tempting to target the nearby star forming regions, which are typically twice as far as the bulk 
of objects usually combed for planets by direct imaging. Probing these interesting reservoirs sets additional constraints that we review 
in this paper by presenting the planet search that we initiated in 2008 around the disk-bearing T Tauri star IM Lup, which is part of 
the Lupus star forming region (140-190 pc). 

Methods. We show and discuss why age determination, the choice of evolutionary model for both the central star and the planet, 
precise knowledge of the host star proper motion, relative or absolute (between different instruments) astrometric accuracy (including 
plate scale calibration), and patience are the key ingredients for exoplanet searches around more distant young stars. 
Results. Unfortunately, most of the time, precision and perseverance are not paying off: we discovered a candidate companion around 
IM Lup in 2008, which we report here to be an unbound background object. We nevertheless review in details the lessons learned 
from our endeavor, and additionally present the best detection limits ever calculated for IM Lup. We also accessorily report on the 
successful use of innovative data reduction techniques, such as the damped-LOCI and iterative roll subtraction. 



Key words. Planet-disk interactions ~ Stars: variables: T Tauri - 
planetary systems - Techniques: high angular resolution 



(Stars): planetary systems - Stars: individual: IM Lup - Infrared: 



- 1 — \ 

X: 
^3: 



1. Introduction 

Direct imaging constitutes an attractive technique for exoplanet 
detection as it provides straigh tforward means to c haracterize 
planets and their h ost system dAbsil & Mawet 2010) through , 
e.g., orbital motion dSoummer et aLluOl ll IChauvin et al. 20121). 
spectro-photometry of planetary atmospheres (iJanson et al.l 
[2010; Galich er et al 2011; Bonnef oy et al. 201 1]), or planet-disk 
interactions dLagrange et al.il2012l) . Direct imaging has also the 
potential of understanding and bridging the gap between the 
population of extremely close planets discovered by radial veloc- 
ity or transit techniques and the free floating plane t s discovered 
by m icrolensing observations tSumi et al. 1201 it lOuanz et al.l 
I2OI2I) . Indeed, many exoplanet candidates directly imaged so 



* Based on the ESO observing programs 380.C-0910, 084.C-0444, 
287.C-5040; and HST observing program 10177. 
** FNRS postdoctoral researcher. 



far have projected distances up to several hundreds of AU. On 
the other hand, some free floating low-mass objects have been 
found to be kine matically associated at projected distances of 
thousands of AU (ICaballero et al.ll2006h . This raises the ques- 
tions of their formation and the very definition of planets, on 
which direct imaging is key to shed more light. 

However, imaging extra-solar planets around other stars con- 
stitutes a multiple challenge, and the practical hurdles are nu- 
merous. First of all, the angular separation between planets and 
stars is very small (e.g. <500 mas for a 5-AU distance at 10 
pc), usually requiring diffraction limited capabilities on 8-meter 
class telescopes Q. Second, the contrast between a planet and 
its host star ranges from ^ 10"^ for hot giant planets in the in- 



' Note the exception presented in Sera bvn et al.l ( |201(]|) . who showed 
a snapshot of 3 out of the 4 planets of HR 8799 taken with an adaptively- 
corrected 1 .5-meter telescope and a ne xt-generation vector vortex phase 
mask coronagraph jMawet et al.ll201(ih . 
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frared to ^ lO '" for Earth-Uke planets in the visible. The con- 
trast issue requires exquisite image (hence wavefront) quality 
to feed coronagraphic devices, most of the time very special- 
ized observing strategie s (e.g., angular differential imaging or 
ADI. lMarois et aT l2006l) . and corresponding data reduction tech- 
niques such as the Locally O ptimized Combination of Images 
(LOCI.' Lafreniere et al.ll2007n . 

Once a faint point source has been detected around a star, 
pointing to the potential discovery of a companion candidate, 
precise differential astrometric monitoring of the latter needs to 
be carried over a sufficiently long time so that the stellar proper 
motion overcomes the astrometric precision of the detected ob- 
ject by a suffi cient margin (if the object is b ound, it moves with 
its host star). iNeuhauser & Schmidt! (1201 2h also argues that a 
spectrum, when possible (proximity to the host star often pre- 
vents to take clean uncontaminated spectra), can determine the 
spectral type and temperature of the companion, and thus in- 
dicates a planetary mass or sub-stellar body, but still possibly 
a cool background object. Both tests might sometimes be nec- 
essary, especially when targeting young associations where ob- 
jects can potentially share common proper motion, likely to be 
small and rather uncertain (at the distance of star forming re- 
gions), making this astrometric process more difficult and the 
required time baseline longer. The T Tauri star ScoPMS 214 is 
a typical example, where a candidate companion was shown to 
share common proper motion, but w as spectr oscopically iden- 
tified as a foreground M dwarf ( Met chev & Hil lenbrand 2009). 
In young associations, the probability for small and/or shared 
proper motion is thus significant. A third possibility for confirm- 
ing the bound character of the companion is the detection of the 
orbital motion, but that implies that the candidate is on a rea- 
sonably tight orbit (period < 1000 years), in order to be sampled 
with sufficient accuracy over a time baseline of a few years. 

Most of the objects imaged so far are orbiting young 
stars (see exoplanet.eu fo r a tho rough and up-to-date list, 
and INeuhauser & Schmidt! (!2012!) for a recent detailed re- 
view). Youth is the current bias of high contrast imaging, 
as short period, inclination or distances (orbital and/or par- 
allactic) are the biases of radial velocity, transit and micro- 
lensing techniques, respectively. Indeed, the thermal radiation 
of young exoplanets peaks in the near-infrared, making them 
more easily detectable by several orders of magnitude than 
if we were to observe them in reflected light in the visible. 
Since the detected emission comes from the intrinsic ther- 
mal radiation of the planet, its physical properties (tempera- 
ture, mass and radius) can only be inferred based on cool- 
ing track models, which critically depend on age and forma- 
tion mechanis ms/history (Allard et al. 2003; M arlev et al.!!200'7t 
!Fortnev et al.ii2008: Spie gel & Burrows 2012). Due to this very 
high sensitivity to initial conditions, deriving the nature of low- 
mass objects and young planets is indeed more problematic than 
for older ones, especially for long-period companions where 
the dynamical mas s is more difficult to infer. For instance, 
!Stassun et al.l (!2006h found a young eclipsing brown-dwarf bi- 
nary in which the cooler object is the more massive one, which 
is very surprising, because most theoretical models predict that 
a brown dwarf of a given mass will at all times be warmer than 
a lower-mass brown dwarf of the same age. lStassun et al.l(l2008h 
also found a binary where both stars have the same mass within 
2% but their surface temperatures differ by 300K. Therefore, dy- 
namical mass, hence astrometric precision and proper motion 
knowledge, age, together with distance, become key parame- 
ters that will determine the final precision and confidence on the 
companion physical characteristics. 



Table 1. Fundamental properties of IM Lup and associated opti- 
cally thick circumstellar disk. 



Properties 


Value 


Names 


IM Lup, Sz82, PDS75 




IRAS 15528 - 3747 


Spectral type 


MO 


Temperature 


3900K 


Class 


CTTS/WTTS 


Age 


0.5 - 1.75 Myr" 


Association 


Lupus 


Distance 


140- 190 pc" 


V mag 


11.96 


H mag 


8.089 


K mag 


7.739 


Lp mag 


7.29 


Disk radius*" 


0.32- 400 AU 


Disk inclination'' 


^ 50° 


Disk-to-star mass ratio'' 


^ 0.1 


Grain sizes'" 


0.03 - 3000 



See the discussion in Sect. [7] 

'' Taken from the best fitted model presented in lPinte et al.l ( 1200 8h . 
- References for all other values are in the text. 



In the present paper, we discuss these issues in details, and il- 
lustrate them with our search for planetary candidates around the 
young T Tauri star IM Lup. The paper is organized as follows: 
after presenting the IM Lup system in Sect. |2] we describe our 
VLT/NACO 2008 observations and discovery of a putative can- 
didate companion in Sect. [3] Then, in Sect. H) we present our re- 
analysis of the pre-discovery data taken with HST/NICMOS in 
2005, followed by epoch 3 and 4 images taken with VLT/NACO 
in 2010 and 201 1 in Sect.|5] In Sect.|6] we discuss our discovery 
in details and underline the difficulties of planet searches around 
this type of objects, in terms of astrometry and related proper 
motion analysis. We also establish benchmark detection limits 
around IM Lup after reestablishing its age (Sect. |7]i, before con- 
cluding in Sect. |8] 



2. IM Lup: a young T Tauri star with a massive 
circumstellar disk 

IM Lup (TableD is a young MO (T= 3900 K) T Tauri star (TTS) 
with an equivalent width of the Ha emission known to vary 
from 7.5 to 21. 5A, confirming its status as a borderline weak- 
line/classical TTS. Part of the Lupus association (140-190 pc, 
see Sect. |7]i, it is one of four young stellar objects in the small 
'^CO(l - 0) Lupus 2 core near the extreme T Tauri star RU Lup 
(Tachihara et al. 1996). Our age estimate described in Sect. |7] 
yields 0.5-1.75 Myr. 

Despite the low accretion-related act ivity of IM Lup 
(!Reipurthet al.1! 19961: !Wichmann et alJ!l999!), long wavelengfli 
observations from the millimeter (^Nuernbe rger et al.l !l997!: 
van Kempen et al. 2007; Lommenetal. 2007) to the infrared 
(iPadgett et al.l 120061) reveal ample evidence for gas-rich cir- 
cumstellar material in the system. IM Lup's protoplanetary 
disk scattered light was imaged in 1999 in the visible with 
HSTAVFPC2 (PI: Stapelfeldt, Prog. ID 7387). It was followed 
in the near-infrared by HST/NICMOS images obtained in 2005 
(PI: G. Schneider, Prog. ID 10177). An extensive modelin g study 
of the IM Lup disk was performed bv !Pinte et al.! (!2008h . using 
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Table 2. Observing log for IM Lup and the reference stars used in this work. 



Target 


Prog.ID 


a (J2000) 


6 (J2000) 


Filter 


UT date 


Exp. time 


Tel./instr. 


Strat.^ 


Strehl (%)" 


IMLup 


10177 


15h 56' 09" 


-37° 56' 06" 


F160W 


29/03/05 


1350i- 


HST/NICMOS 


ADI 




IMLup 


380.C-0910(A) 


15h 56' 09" 


-37° 56' 06" 


Ks 


29/03/08 


1350i- 


VLT/NACO 


RDI 


58 + 5 


CD-37 8989 


380.C-0910(A) 


13h 54' 27" 


-38° 14' 54" 


Ks 


29/03/08 


1200i 


VLT/NACO 




50 ±3 


CD-35 9033 


380.C-0910(A) 


13h 50' 00" 


-36° 33' 40" 


Ks 


30/03/08 


1200i 


VLT/NACO 




48 ±3 


IMLup 


084.C-0444D 


15h 56' 09" 


-37° 56' 06" 


Lp 


19/04/10 


1350J 


VLT/NACO 


RDI 


80 ±5 


LHS3286 


084.C-0444D 


17h 23' 49" 


-32° 15' 16" 


Lp 


19/04/10 


13505 


VLT/NACO 




85 ±5 


IMLup 


287.C-5040(A) 


15h 56' 09" 


-37° 56' 06" 


Ks 


22/07/11 


13505 


VLT/NACO 


I 


40 ±3 


IMLup 


287.C-5040(A) 


15h 56' 09" 


-37° 56' 06" 


H 


25/07/11 


13505 


VLT/NACO 


I 


30 ±3 



" Observing strategy: ADI (angular differential imaging), RDI (reference star differential imaging), I (imaging). 
The Strehl ratio was measured on the reduced image, or on the acquisition PSF (for saturated or coronagraphic data). Seeing and coherence 
time (to) conditions in the Visible (^ 0.5 yum) were very good for all observations, typically 0'/6 - 0'.'8, and 3-8 ms, respectively. 

- Note that all targets were observed at airmass 1.1, except for HST of course. 



multi-wavelength spectro-photometry and images in a global fit 
with a 3D radiative transfer model, and led to quantitative evi- 
dence for dust processing an d evolution i n the d isk. 

A more recent paper by iPanic et al] (12009') presents SMA 
(Submillimeter Array) observations showing a break in the gas 
and dust surface density of the IM Lup disk, seen to exten d much 
furth er than the 400 AU outer edge determined by Pint e et al] 
(1200 8). One of the proposed explanations for the break is a com- 
panion body near the break at 400 AU. Indeed, a companion 
of 1 Mjuo could open a gap in the disk and affect its spread- 
ing. Panic et an (l2009h however argue that no candidate compan- 
ion_at this separation is visible in the HST image of .Pinte et aL 
(l20(j8h . 

3. Discovery of a candidate companion with 
VLT/NACO in 2008 

As part of a coronagraphic study of young stars (prog. ID 380. C- 
0910(A), PI: Mawet), we observed IM Lup in March 2008 with 
NAOS-CONICA, the adaptive optics (Nasmyth Adaptive Optics 
System) and near-infrared spectrograph and imager of the Very 
Large Telescope (VLT). 

3.1. Observing strategy for the 2008 discovery data set 

For our discovery image in 2008 , we used the fou r-quadrant 
phase-mask (FQPM) coronagraph (iRouan et al ] l2007h in the Ks 
band. The FQPM is a phase-mask coronagraph applying a n 
phase shift between adjacent quadrants. The starlight, when cen- 
tered on the FQPM cross-hair, undergoes a destructive interfer- 
ence upon propagation to the downstream pupil plane. 

All of our frames for 2008 were taken with the Ks filter and 
the S13 camera (13.27 mas/pixel). This fine sampling (4 pixels 
per resolution element A/d, where A is the observing wavelength 
and d the telescope diameter) was necessary to center the target 
star on the FQPM cross-hair precisely. The main calibrator stars 
were carefully selected to present roughly the same V-K color 
as the target (Table |2]i. Matching the V magnitudes is important 
to ensure similar AO corrections between the target and refer- 
ence stars, as the visible wavefront sensor of NAOS is mostly 
sensitive at V and R. Ks magnitudes have also to be matched 
to ensure SNR matching for the quasi-static speckles. Also, to 
avoid flexure-induced semi-static speckle variations as much as 
possible and to ensure a consistent telescope orientation with re- 
spect to the instrument between the target and the reference, the 
calibration stars were chosen and observed at the same paral- 



lactic angle as the target star This condition was met on a best 
effort basis since the availability of a suitable reference fulfill- 
ing the set of constraints is never guaranteed, which is one of 
the drawbacks of the reference star differential imaging strategy 
(RDI). 

The observing conditions for IM Lup and the reference stars 
were very good with a visible seeing between 0'.'6 and 0'/8. The 
total integration time was about 1350 sec and 1200 sec for each 
target (see Table|2|. Respecting the consideration discussed here 
above to calibrate time-dependent PSF variations (speckle), we 
acquired coronagraphic images of the reference star 90 min af- 
ter the science target at roughly the same parallactic angle. To 
reduce drift and pupil rotation, the centering was checked and 
corrected every 80 sec. 

3.2. Data anaiysis for tiie 2008 discovery data set 

The data reduction proceeds as follows. NACO coronagraphic 
acquisition template moves the telescope alternatively between 
a fixed object position and a jittered set of sky positions which 
are median combined and subtracted to the object, removing the 
background and dark contributions at the same time. The nor- 
malization of the resulting image with the flat provides the first 
stage of data processing. However, this treatment does not com- 
pletely remove the electronic noise of the readout process nor 
the two electronic ghosts which are characteristic of CONICA. 

The second stage of the data reduction process consists in co- 
adding of the images with a sub-pixel centering procedure. For 
that, we applied a hybrid method which correlates the centroid 
of the unsaturated coronagraphic pattern with a pre-computed 
look-up table based on a diffraction model of the FQPM. Using 
this sophisticated method, we routinely achieve a centering pre- 
cision of cr = 0.1 pixel or 1.4 mas rms. Despite the coronagraph 
starlight 10-fold attenuation, scattered starlight still dominates 
the extended source flux. Since the primary objective of our orig- 
inal program was disk imaging, and given the relative novelty of 
the ADI technique at the time of the first observations, and the 
fact that ADI is not an optimal strategy for nearly face-on disk 
(Milli et al. 2012, accepted to A&A), we chose to use classical 
RDI. 

The final VLT/NACO 2008 image (see Fig. [T] top middle) 
was obtained using an enhanced version of the LOCI algorithm 
("Lafreniere et al.ll2007h . In short, LOCI finds the optimal linear 
combination of reference frames (here from the two reference 
stars CD-37 8989 and CD-35 9033) to minimize the noise in a 
given zone of the target image. The process is repeated until the 



3 



D. Mawet et al.: Direct imaging of extra-solar planets in star forming regions 



NACO 2008 (Ks) - RDI 



NACO 2008 (Ks) - RDI - dLOCI 



NICMOS 2005 



'J 

k«4 



SW arc 



SW arc 



1" 



NACO 2010 (Lp) 



1" 




Fig. 1. The upper strip presents the simple reference star subtracted Ks-band image (left) and the damped-LOCI processed image 
(middle) obtained with NACO in 2008, which show a very faint off-axis point source {iuks - 19.1) to the northeast of IM Lup. The 
upper right image presents the 2005 F160W NICMOS observations, using an iterative roll-subtraction technique. The lower strip 
presents the NACO Lp-band image of 2010 (left), the 201 1 NACO H-band image (middle), and the 201 1 NACO Ks-band image 
(right). The off-axis point source is detected in all images (identified with an arrow). 



area of interest in the target image is completely reduced. LOCI 
in its original form was conceived to find point sources, and has a 
known tendency to attenuate signal from extended sources such 
as circumstellar disks. However, this defect of the generic LOCI 
algorithm can be brought under control by a fine tuning of the 
geometrical parameters such as the size of the optimization zone, 
the number of reference frames used in the correlation matrix, 
as well as the introduction of a damping parameter (Lagrange 
mult ipher) to balance fl ux conservation with noise attenuation 
as in lPuevo et all (|2012|) . 



3.3. 



Candidate companion and new image of tlie disk 



In our 2008 Ks-band dataset (Table |2]l, we detected a very faint 
companion at a signal-to-noise ratio ( SNR) of ^ 1 (Fig.[T] top 
left). Applying the damped-LOCI of iPuevo et alJ (1201 2 ). both 
the companion and the disk SW arc are detected simultaneously 
with minimum contamination from starlight scattering (Fig. [1] 
top middle). For the first time, we note that damped-LOCI, orig- 
inally invented to detect point-source in multi-spectral data, can 
be successfully applied to the detect ion of circumste l lar dis ks, 
improving upon the original LOCI of iLafreniere et al.l (l2007h . 



We performed relative astrometry on the final reduced im- 
age, using centroiding and a specific pre-computed look-up table 
for the star position behind the FQPM (see Sect. 13.2b . and gaus- 
sian fitting for the candidate companion (with a subpixel preci- 
sion). The companion is located to the North-East of IM Lup, 
at a radius of ^ 1'.'8, and a position angle (PA) of ^ 58° (see 
Table©. 



We then performed aperture photometry using the function 
APER of IDL, and found a relative Ks magnitude of IQ.l^g j^. 
Naively assuming that the point source is physically associated 
to IM Lup (Sect. |6]), this corresponds to an absolute Ks magni- 
tude of 13.0^g g, where the uncertainty is mostly due to the poor 
knowledge of the star distance (140-190 pc, see Sect.|7]). 



Note that the disk is detected as an arc to the SW, extending 
up to ^ r.'8 along the major axis, and ^ 1'.'3 along the minor 
axis, consistent with previous HST observations. We also mea- 
sure a position angle of 140° + 10° f or the major axis , which is 
consistent with the value reported in lPinte et al.l (l2008h . 
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Table 3. Astrometry of the point source to the NE, along with its relative and absolute photometry. 



Data set 


Filter 


SNR 


Aa (") 


AS (") 


m 


M" 


NICMOS 2005 
NACO 2008 
NACO 2010 
NACO 2011 
NACO 2011'' 


F160W 
Ks 
Lp 
Ks 
H 


3 

10 

^ 2 
^ 10 
^ 5 


1'.'526 + a'043 
1'.'541 ±0'.'007 
1'.'570 ± a'020 
1'.'575 + a.'009 
1'.'594 ± a.'025 


a'937 + a.'043 
a'936 ± a.'007 
a.'978 ± a.'020 
1'.'002 + a.'009 
a'996 ± a.'025 


19 3+"" 

19 r"- 

17.8!«J 
19 25+"-^ 


13.0!°-^ 
11. 7+'-* 

13 r^-^ 

1 ^ 9+0.6 

^■^•^-04 



" Assuming physical association. The error bars include the distance uncertainty (140-190 pc). 

'' The error bars are large due to the presence of the diffraction pattern of the telescope spiders close to the point source. 



4. Pre-imaging of the candidate companion with 
HST/NICIVIOS in 2005 

Following the NACO discovery, we retrieved the archival 
HST/NICMOS data of IM Lup obtained in 2005 in the F160W 
filter (PI: G. Schneider, Prog. ID 10177). The reason for the non- 
detection of the co mpanion by successive groups and notably by 
iPinte et al.l ( |2008|) is rather simple: they primarily aimed at char- 
acterizing the extended structure, and applied corresponding soft 
data reduction techniques (single reference star subtraction). 

Taking advantage of the availability of two images taken at 2 
different roll angles separated by 30 degrees , we used the itera - 
tive roll-subtraction algorithm introduced in iHeap et aP (l2000h . 
and described in Krist et al.l 02010 ). A simple roll-subtracted im- 
age contains a positive image of a companion and a negative one 
rotated by an angle equal to the difference between the telescope 
orientations. Ideally, these two would be combined to form a sin- 
gle, positive image. Instead of directly subtracting images from 
the two orientations and then trying to combine the results, an 
iterative technique can be used that solves for those portions of 
the two unsubtracted images that appear static on the detector 
(i.e., the PSF) and those that appear to rotate as the telescope 
rolls (i.e., the sky, including any companion or disk). 

This method, which has been successfully used for /J Pictoris 
jHeap et al. .2000) and HD 207129 ( Kiist et al. 2010), clearly 
(^ 3cr) reveals the companion in the NICMOS data as well 
(Fig. [H top right). The position of the star behind the opaque 
Lyot coronagraph of NICMOS was determined by two different 
methods. The first position was obtained by iteratively shifting a 
Tiny Tim model (Ki'ist.l993) PSF by subpixel amounts until the 
residuals were visibly minimized. The second position estimate 
was obtained by measuring the intersection of the diffraction 
spikes of the secondary mirror support structures. Both meth- 
ods yield similar result within half a pixel, which we considered 
as a systematic error folded in the error bar calculation (TableO. 

5. Epoch 3 (2010) and 4 (2011) with VLT/NACO 

Given the insufficient time baseline and astrometric precision 
of the NICMOS data point, we waited for a few years and re- 
observed IM Lup in 2010 (Lp band) and in 2011 (H and Ks 
band), using NACO again (Table |2|i. The Lp-band data only 
marginally (SNR^ 2) shows the companion (Fig.[T] bottom left), 
with an estimated Lp magnitude of 17. 8;^^ 5. On the other hand, it 
is easily detected in H and Ks bands in the 201 1 data (Fig.[T] bot- 
tom middle and right). The strategy we chose for the most recent 
data set was to perform simple saturated-unsaturated imaging in 
order to enable precise astrometric and photometric analysis. 

Results of the astrometry of the candidate companion rela- 
tive to the host star are presented in Table[3] Note that the relative 
astrometry is somewhat different between both filters. Slight dif- 



Table 4. Proper motions for IM Lup found in major astrometric 

catalogues 



Origin" 






H 




(mas/yr) 


(mas/yr) 


(mas/yr) 


Hipparcos 


-56.66±15.41 


-49.97+7.31 


76±17 


Leeuwen 


-35.5±21.81 


-22.93±14.98 


42±26 


PMS* 


-3±2 


-21 ±2 


21±3 


SPM4 


-12.7+3.9 


-21.5±4 


25±6 


PPMXL 


-2.03±4 


-21.53±4 


22±6 


UCAC3 


-15.4±4.7 


-22.6±5.1 


27±7 



" Column (1) gives the origin of the measurem e nt: L eeuwen: 
re-re duction of Hipparcos data Ivan LeeuwenI |2007|) . PMS 
(Duc ourantetal . 2005), SPM4: the Yale/San Juan SPM4 Catalog 
(2009 http://www.astro.yale.edu/astrom/spm4cat/spm4.html ), 
PPMXL (Roeser et al. 2010), UCAC3 (Zacharias et al. 2010). 

*" Column (2) and (3) give proper motions together with uncertainties. 
Column (4) gives the modulus of proper motion. 

* Adopted in this work. 



ferences are expected, due for instance to differential aberrations 
between filters, and the difference in data quality (the Strehl ra- 
tio is naturally lower in the H band). However, it appears that 
the H-band astrometry is affected by the presence of the diffrac- 
tion pattern of the telescope spiders (secondary mirror support 
structure), close to the point source. For these reasons, we only 
retained the Ks band astrometry in our final proper motion anal- 
ysis of Sect. 16. II which also corresponds to the filter of the plate 
scale calibration described below. 



6. Discussion 

In this section, we elaborate on the difficulty of exoplanet candi- 
date confirmation and characterization for young distant stars. 

6. 1 . Astrometry and proper motion 

To test the hypothesis of an object linked to IM Lup we can use 
astrometry. If the companion is linked to IM Lup, no significant 
variation of their angular separation should be observed along 
time due to proper motion, since they will be co-moving. The 
only cause of variation would result from the orbital motion of 
the companion around the center of mass of the system. The ap- 
parent angular separation of ^ I'.'S corresponds to a deprojected 
physical separation of about 350-480 AU at 140-190 pc. We can 
expect a planet orbiting at such a distance from its central star 
to have extremely long periods (several thousand years at the 
minimum) so that its orbital motion is not detectable in our as- 
trometric observations spread over 7 years. 
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6.1 .1 . Prior astrometric calibration of CONICA 

To calibrate the NACO plate scale and detector orientation in a 
consistent and precise way between the 2008 and 2011 epochs, 
we used the star clusters Theta Orionis/'Trapezium" (2008) and 
Omega Centauri (2011). The reference positions of the stars in 
the each cluster were derived in a different program (Montagnier 
et al. 2012, in preparation). For each cluster, many images were 
taken at various positions and orientations to establish a distor- 
tio n solut ion of the N ACO plate-scale (see the method described 
by [Anderson & Kinj ( l20()3ft ). The linear terms of the distortion 
(detector axes orientations and pixel dimensions) were then de- 
rived by observing appulses of transneptunian objects. Finally, 
the star positions were derived with an accuracy better than 1 
mas (about 200 object for the Omega Centauri cluster field, and 
about 50 for the Theta Orionis cluster). On the calibration im- 
ages of the 2 epochs needed in our astrometric analysis, the po- 
sition of the centroid (s^, Sy) of each non-saturated star (about 
20 stars in the Theta Orionis 's 2008 epoch, and about 40 for 
the Omega Centauri's 201 1 epoch) on the reduced image of the 
field is measured in pixels; these values are then compared to 
the position on the sky (Pa,Pv) in arcseconds with the following 
equations: 

PjC ^ X0+ P_,Sj, cos 0, - PySy sin Oy (1) 

Pv = yo + PxSx sin Oj, + pySy COS Oy (2) 

where (/?,,) is the plate scale along the x- (y-) axis, 9x (Oy) 
are the orientations of the detector on sky along the x- (y-) axis, 
and xo and yo are offsets giving the correspondence between 
the absolute positions (it is only used to solve the equation). A 
Levenberg-Marquardt minimization is then applied to find the 
plate scale solution. Using this calibration method, the final pre- 
cision on the plate scale is ^ 50/ias, and +0.05° on the detector 
orientation. Note that this accuracy is only one item in the er- 
ror budget of the final astrometric precision, which depends on 
many other terms all summed quadratically, such as the preci- 
sion of the star position determination, fit of the off-axis location 
(itself dependent on the SNR), etc. 

6.1.2. Proper motion of IM Lup 

We have considered the proper motions of IM Lup available in 
the literature (see Table lU to test our hypothesis of co-moving 
objects. We notice that the measured proper motions vary in a 
large range from one author to the other The p utative binary na- 
ture of IM Lup wro ngly reported by Hippa rcos dWichmann et al] 
ll998HLasker etal . 1996; Kharchenko & R oeseijl2009h indicates 
a potential disturbance in the measurements that led to a poor 
astrometry. Indeed, our observations do not reveal the second 
component down to a magnitude Mks - 19, indicating that 
the Hipparcos detection, and other reports of the binarity might 
have been potentially contaminated by the presence of the opti- 
cally thick circumstellar disk of IM Lup, as already suggested in 
iPinte et al.1 (120081) . 

Other authors have measured the proper motion of this un- 
resolved object with unequal precisions, leading to a consensual 
value in declination (~22.4 mas/yr) while the proper motion in 
right ascension varies largely between authors (from -3 to -15.5 
mas/yr for values with reasonable precisions). The origin of such 
discrepancies is difficult to pinpoint since proper motion quality 
is not only related to the time base but also to the number of dif- 
ferent epochs of observation and evidently to the quality of each 
epoch measurement. In the case of pre-main sequence stars the 
situation is even more difficult since depending on the target, the 



object may be embedded into a dust and/or gas cloud (which is 
the case for IM Lup) perturbing the photo-center measurement. 
Moreover the morphology of the cloud may vary with time and 
lead to variable photo-centers at different epochs. In the data pre- 
sented in TablelH we chose to adopt the third one (PMS), but the 
four last values may be considered (PMS, SPM4, PPMXL and 
UCAC3) for the astrometric test of co-moving objects. 

6.1 .3. Bound or not bound? 

With these proper motions of IM Lup, we would expect a 
background source to have moved by 61 +9 mas with respect 
to IM Lup between the two observations (NICMOS-2005 and 
NACO-2008). Such a motion is not detected within our error 
bars, meaning to first order that the companion is likely co- 
moving with IM Lup (Fig. |2] left). However, the average SNR 
on the NICMOS detection, the very slow proper motion and 
the large astrometric uncertainty mentioned above, do not al- 
low us to firmly and definitely conclude on its bound charac- 
ter. Note that a galactic star count model for the direction toward 
IM Lup ( Girardi et al.l2005h yields a surface density of stars with 
19 < Ks < 19.5 of ^ 2.2 x 10^ per square degree. This makes 
the chance of a random background source being located within 
l'.'8ofIMLup:^ 17%. 

Summer 201 1 was the first opportunity to firmly get closure 
on the bound aspect of this discovered candidate. We then used 
the epoch 4 NACO observation to redo the common proper mo- 
tion astrometric analysis. This time, since the analysis is based 
on a single, well calibrated instrument, our astrometric precision 
can be trusted down to a conservative ^ 10 mas per coordinate. 
With a time baseline of 1210 days, the background object should 
have moved by 68 + 10 mas with respect to IM Lup, which is 
about the observed variation of separation (74 + 20 mas) in the 
same direction. We conclude that the candidate companion is 
likely to be a background object, and is therefore not associated 
with IM Lup (Fig. El right). 

6.2. Probable nature of the point source 

Based on the combined H, Ks, Lp photometry, we analyzed the 
SED of the likely background object to verify that it is consis- 
tent with a blackbody. For that, we first checked that the extinc- 
tion in the direction of IM Lup (^ 25° from the galactic center 
bulge) is very small and can indeed be neglect ed in the near- 
infrared: Ah - 0.3, Ak ^ 0.2, an d Al„ ^ 0.1 jSchlegel et al.1 
ll998l:ISchlafly & Finkbeineil201 1 ). The SED would be compat- 
ible with many possible stellar objects. For instance, the fit to a 
3000 K blackbody is satisfying, with the Lp-band point falling 
only a Uttle more than one sigma above the model. Any black- 
body warmer than about 3000 K would actually fit the SED in 
almost the same way. The reduced of the fit is ^ 2 when both 
the size and temperature of the blackbody are simultaneously fit- 
ted. If the background object is actually an M5V star, it would 
be located far away in our galaxy (^ 4kpc). 

Finally, confirming the stellar nature of the point source, we 
verified that the H-Ks color derived from our data is not con- 
sisten t with usual evolution ary models for p lanetary-mass ob- 
jects dBaraffeet aL 2003; Fo rtnev et al.ll2008l) . which are gener- 
ally much redder. Note that the Ks-Lp color of ^ 1.5 - 2 we 
measured is marginally consistent with both hot start and core 
accretion models, but the error bars are such that the Lp photo- 
metric point is not significant, hence constraining. 
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Fig. 2. Left: astrometric analysis of the candidate companion proper motion between 2005 and 2008, using NICMOS and NACO. 
Right: astrometric analysis of the candidate companion proper motion between 2008 and 201 1 (and the marginal Lp-band detection 
of 2010), using NACO only. The continuous line represents the combined parallactic motion and the proper motion of the back- 
ground object in the reference frame of the primary target. The filled dots represent the positions of the candidate companion with 
error bars at the different epochs while the empty dots ("IF BKGD") represent the positions of a putative background object at the 
same epochs. 



7. Age, distance, and detection limits 

Preliminary age estimates f or IM Lup range from 0.1 Myr to 
10 Myr (Hu ghes et al.lll994l) . To reduce the uncertainty associ- 
ated with this large range, we re-estimated the age of IM Lup 
as follows. We placed the object on a Hertzsprung-Russell dia- 
gram (HRD). The effective tem perature w as given by the spec- 
tral type (conversion from L uhman et alJ ([2003.) for M dwarfs). 
I and J magnitudes, which are not too much affected by ac- 
cretion nor disk emission, were converted to bolometric mag- 
nitudes based on bolometric corrections and intrinsic colors of 
iKenvon & Hartmanni(f 19951) . 

The bolometric luminosity was then deduced using an 
estimated distance most probabl y comprised betw een 140 pc 
(iHughes etalJlT994l) and 190pc Iwichmann e taT 19980, and 
corrected for extinction u sing an Av = 0.5 (IPinte et al.. .2008.) 
with the law presen ted in Draina (l2003h . W e then used evolu- 
tionary models from'Baraffe et al.' (^1998) and' Siess etaP (l2000l) 
to draw isochrones and evolutionary tracks in the HRDs and to 
interpolate for the observed object. The age estimation was per- 
formed independently for I and J photometry, and then folded 
into error bars. We arrived at the following estimates: for a dis- 
tance of 140 pc, we get an age range of 0.8-1.75 Myr, while for 
a distance of 190 pc, we get 0.5-1 Myr 

Our 2008 FQPM data set taken in the Ks band features the 
best contrast ever achieved around IM Lup. This data set is there- 
fore suitable to derive detection limits for IM Lup (Fig. [3]|. We 
proceeded as follows. For increasing angular separations, we de- 
rived the standard deviation in annuli 1 resolution element wide. 



^ Recently confirmed by Galli et al. 2012 (in preparation), who mea- 
sured a kinematic distance of 179 pc for IM Lup. 



This profile was then multiplied by 5 to derive the 5cr detection 
limit associated with the corresponding data set. 

In Fig. [3] we also overplot the level of contrast for two planet 
masses: 1, 2 Mjup, assuming a me dian age of ^ 1 Myr and con- 
sidering the COND03 model from lBaraffe et al. ' (2003). Thanks 
to the young age of IM Lup, our (model-dependent) detectability 
limits are excellent, down to less than one Jupiter mass beyond 
one arcsecond, and less than two Jupiter masses beyond 0'.'2. 
Note that the "core accretion" model of Fortney et al. (2008|) 
yields much higher masses, reflecting the large uncertainties 
still p laguing evolutionary models for early ages dMarlev et al.l 
120071) . 



8. Conclusion 

This paper presented a planet search we conducted with 
VLT/NACO around the young T Tauri star IM Lup between 
2008 and 2011, using a pre-discovery image obtained with 
HST/NICMOS in 2005. IM Lup is the perfect prototype system 
for planet search since it has a massive optically thick circum- 
stellar disk, likely at the stability limit. It also features a break in 
the gas and dust density at about 400 AU, which could indicate 
the presence of a Jupiter-mass body at the location of the discon- 
tinuity. A candidate companion was detected by NACO in 2008, 
and also seen in the 2005 HST/NICMOS data. 

The candidate companion is located to the North-East of 
IM Lup, at a radius of =^ V.'S, and a position angle (PA) of ^ 58°. 
Tentatively and naively assuming association, this corresponds 
to a de-projected physical separation of about 350-480 AU at 
140-190 pc. With our redetermined age of about 1 Myr, the mass 
of the putative off-axis companion using the usual "hot start" 
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Fig. 3. 5cr detection limit around IM Lup, derived from the Ks- 
band 2008 data, assuming an age of ^ 1 Myr and the COND03 
model from Baraffe et al. (2003). The dashed curves show the 
limits for distances of 140 and 190 pc, respectively. The dashed 
lines shows the level of contrast for two planet masses: 1, 2 Mjup. 
The black circle and associated error bar shows the point-source 
Ks band absolute magnitude if associated. 

evolutionary models (iBaraffe et alJ 120031: iFortnev et aP l2008l) 
would be between 1 - 2Mjup. 

However, and unfortunately, the candidate was later on 
proven to be a background object based on the NACO 201 1 ob- 
servations, a common proper motion analysis and a careful cal- 
ibration of the NACO plate scale and detector orientation. This 
cautionary tale taught us the difficulty of planet search around 
young, distant and obscured stars, where proper motion might 
not be very well constrained, and where the age and distance 
determinations are tricky. 
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